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In  most  microgravity  studies  of  flame  spread,  the  flame  is  assumed  to  be  two-dimensional,  and  two- 
dimensional  models  are  used  to  aid  data  interpretation.  However,  since  limited  space  is  available  in  mi¬ 
crogravity  facilities,  the  flames  are  limited  in  size.  It  is  important,  therefore,  to  investigate  the  significance 
of  three-dimensional  effects.  Three-dimensional  and  two-dimensional  simulations  of  ignition  and  subse¬ 
quent  transition  to  flame  spread  were  performed  on  a  thermally  thin  cellulosic  sample.  Ignition  occurred 
by  applying  a  radiant  heat  flux  in  a  strip  across  the  center  of  the  sample.  The  sample  was  bounded  by  an 
inert  sample  holder.  Heat  loss  effects  at  the  interface  of  the  sample  and  the  sample  holder  were  tested  by 
varying  the  thermal-physical  properties  of  the  sample  holder.  Simulations  were  also  conducted  with  samples 
of  different  widths  and  with  different  ambient  wind  speeds  (i.e.,  different  levels  of  oxygen  supply). 

The  width  of  the  sample  affected  both  the  duration  of  the  flame  transition  period  and  the  post-transition 
flame  spread  rate.  Finite  width  effects  were  most  significant  when  the  ambient  wind  was  relatively  small 
(limited  oxygen  supply).  In  such  environments,  the  velocity  due  to  thermal  expansion  reduced  the  net 
inflow  of  oxygen  enough  to  significantly  affect  flame  behavior.  For  a  given  sample  width,  the  influence  of 
thermal  expansion  on  the  net  incoming  oxygen  supply  decreased  as  the  ambient  wind  speed  increased. 
Thus,  both  the  transition  and  flame  spread  behavior  of  the  three-dimensional  flame  (along  the  centerline) 
tended  to  that  of  the  two-dimensional  flame  with  increasing  ambient  wind  speed.  Heat  losses  to  the  sample 
holder  were  found  to  affect  the  flame  spread  rate  in  the  case  of  the  narrowest  sample  with  the  slowest 
ambient  wind. 


Introduction 

Many  flame  spread  studies  over  a  solid  fuel  surface 
in  microgravity  have  been  conducted  with  the  mo¬ 
tivation  of  improving  our  understanding  of  the  phys¬ 
ics  relevant  to  fire  safety  in  a  spacecraft.  Almost  all 
these  studies  were  intended  to  observe  the  flame 
spread  behavior  in  a  two-dimensional  configuration, 
and  the  results  were  compared  with  calculated  data 
based  on  two-dimensional  models.  However,  by  ne¬ 
cessity,  these  experiments  were  generally  conducted 
in  spatially  limited  experimental  volumes  which  re¬ 
sulted  in  relatively  narrow  samples  despite  the  origi¬ 
nal  intention  of  a  two-dimensional  experimental  con¬ 
figuration.  For  example,  a  majority  of  experimental 
studies  were  conducted  with  5  cm  width  samples  [  I  — 
3],  but  others  were  with  3  cm  width  [4]  and  6.35  mm 
width  samples  [5],  Future  microgravity  combustion 
experiments  in  space  will  be  carried  out  mainly  on 
the  international  space  station  (ISS),  either  in  a  ded¬ 
icated  combustion  integrated  rack  (CIR)  with  spe¬ 
cific  inserts  for  a  given  experiment  or  in  the  micro¬ 
gravity  science  glovebox.  Drop  towers,  sounding 


rockets,  and  airplanes  will  continue  to  be  used  for 
ground-based  experiments.  The  available  test  cham¬ 
ber  volume  of  these  facilities  will  be  limited,  and 
future  flame  spread  experiments  will  be  restricted  to 
using  relatively  narrow  samples.  It  is,  therefore,  crit¬ 
ically  important  to  determine  whether  these  sample 
widths  will  produce  truly  two-dimensional  flame  be¬ 
havior.  Three-dimensional  effects  on  flame  spread 
along  a  sample  which  is  too  narrow  can,  for  example, 
arise  from  enhanced  oxygen  supply  from  the  outer 
free  stream  and  from  heat  loss  to  a  sample  holder. 
In  this  paper,  transient  two-dimensional  and  three- 
dimensional  simulation  codes  [6]  were  used  to  in¬ 
vestigate  the  effects  of  finite  sample  width  on  flame 
behavior.  These  codes  simulate  the  development  of 
flame  spread  from  ignition,  through  transition  to  es¬ 
tablished  flame  spread. 


Theoretical  Model 

Since  the  theoretical  model  used  for  the  con¬ 
densed  phase  sample  has  been  described  in  Ref.  [7], 
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FIRE — Ignition  and  Flame  Spread 


plane  of  2D  simulation 


Fig.  1.  Schematic  showing  the  cellulosic  sample  and 
computational  domains  of  the  two-  and  three-dimensional 
simulations.  The  two-dimensional  simulation  domain  was 
the  y  =  0  plane.  A  characteristic  char  pattern  from  the 
radiative  ignition  source  is  also  shown.  A  sample  holder 
abutted  the  cellulosic  sample  on  the  y  boundaries.  By 
changing  the  thermal-physical  properties  of  the  sample 
holder,  the  influence  of  heat  loss  to  the  sample  holder  on 
flame  behavior  was  tested. 


only  a  summary  wall  be  given  here.  A  global  ap¬ 
proach  is  used  to  define  both  the  gas-phase  oxidation 
reaction  and  the  thermal  degradation  reactions  of 
the  condensed  phase.  Although  more  detailed  ki¬ 
netic  modeling  is  possible,  the  numerical  solution  of 
the  resulting  system  of  equations  would  be  prohib¬ 
itively  expensive.  The  thermal  degradation  of  the 
cellulosic  sample  is  modeled  by  three  exothermic 
global  reactions  governing  pyrolysis  (slightly  exo¬ 
thermic),  thermal  oxidation,  and  char  oxidation.  The 
evolution  of  the  thermal  degradation  process  is  gov¬ 
erned  by  conservation  equations  for  solid  mass,  cel¬ 
lulose,  char,  and  energy  with  Arrhenius-type  reac¬ 
tion  rates  [8].  The  global  kinetic  constants  and  heats 
of  reaction  were  obtained  by  thermogravimetric 
analysis  and  are  identical  to  those  used  in  Refs.  [6,9], 
As  reported  in  Ref.  [3],  the  samples  thickness  is  0.13 
mm,  area  density  is  5.7  mg/cm  ,  specific  heat  is  0.96 
J/(g  '  K),  and  the  emissitivity  is  0.6.  The  specific  heat 
increases  linearly  with  temperature  and  is  also  used 
for  the  char  and  cellulose.  The  sample  is  assumed  to 
be  thermally  thin  and  of  uniform  composition  with 
depth. 

Radiative  losses  from  the  surface  are  included.  It 
has  not  been  demonstrated  that  radiative  heat  trans¬ 
fer  has  a  significant  effect  on  flames  over  thermally 
thin  fuels  which  are  not  near  extinction.  This  regime 
is  the  focus  of  our  paper.  For  example,  Altenkirch  et 
al.  [5]  have  found  that  for  thermally  thin  fuels,  heat 
transfer  to  the  fuel  surface  drives  the  flame  spread 
and  is  dominated  by  conduction.  Furthermore,  as 
the  sample  width  increases,  the  shape  factor  be¬ 
tween  the  flame  and  the  fuel  surface  ahead  of  the 
flame  will  reach  its  value  for  a  two-dimensional 
flame,  implying  that  the  radiative  flux  to  the  surface 
will  asymptote  to  its  value  for  a  two-dimensional 


flame.  As  will  be  seen  below,  this  is  consistent  with 
the  trend  of  both  the  spread  rate  and  the  duration 
of  the  transition  period  with  increasing  sample 
width.  Thus,  it  is  unlikely  that  including  thermal  ra¬ 
diation  in  the  gas  phase  will  change  the  findings  of 
the  paper. 

The  gas-phase  model  [8]  consists  of  the  conser¬ 
vation  equations  for  total  mass,  momentum  (full  Na- 
vier-Stokes),  energy,  and  species  mass  fractions  (fuel 
and  oxygen)  formulated  in  the  low  Mach  number 
limit.  A  global  one-step  reaction 

[fuel  gas]  +  v02  [oxygen]  — >  product  (1) 

with  an  Arrhenius  rate  is  used.  The  pre-exponential 
factor  in  the  reaction  rate  is  5.0  X  109  cm3/(g  s),  the 
activation  energy  is  67  kj/mol,  the  heat  of  combus¬ 
tion  is  35  kj/g,  and  the  stoichiometric  constant  is 
vGo  =  3.57.  These  values  are  identical  to  those  used 
in  previous  two-dimensional  [3]  and  three-dimen¬ 
sional  [6,9]  studies.  The  average  molecular  weight  is 
assumed  to  be  constant;  temperature-dependent 
molecular  transport  coefficients  are  used.  All  gases 
have  equal  mass  diffusivity  which  is  not  necessarily 
equal  to  the  thermal  diffusivity.  Ignition  is  initiated 
by  an  external  radiative  flux  on  the  sample  surface. 


Numerical  Model 

Both  two-  and  three-dimensional  simulations  were 
performed.  A  sketch  of  the  cellulosic  sample  and  the 
two-  and  three-dimensional  simulation  domains  are 
shown  in  Fig.  1.  The  domain  for  the  two-dimen¬ 
sional  simulations  consists  of  the  y  =  0  (centerline) 
plane.  The  governing  equations  for  the  condensed 
phase  are  solved  using  a  stiff  ordinary  differential 
equation  solver.  A  finite  difference  control  volume 
approach  on  a  staggered  grid  is  used  to  solve  the  gas- 
phase  equations.  A  time-splitting  algorithm  is  used 
for  the  species  and  temperature  equations  since  the 
characteristic  time  scale  of  the  chemical  reaction  is 
relatively  small  compared  to  that  for  diffusion  and 
convection  (high  Damkohler  number)  [8],  This 
method  involves  first  solving  the  species  and  tem¬ 
perature  equations  without  the  convective  and  dif¬ 
fusive  terms  using  a  stiff  ordinary  differential  equa¬ 
tion  solver  over  small  time  steps.  The  global  time 
stepping  is  performed  using  an  alternating  direction 
implicit  scheme  on  the  convective  and  diffusive 
terms.  The  momentum  equation  is  solved  using  a 
projection  method;  the  solution  to  the  Poisson  equa¬ 
tion  for  the  pressure  perturbation  is  found  with  a 
direct  solver. 

There  were  two  assumed  symmetry  planes  in  the 
three-dimensional  simulations.  One  parallel  to  the 
plane  of  the  sample  and  passing  through  its  middle; 
the  other  was  perpendicular  to  the  plane  of  the  sam¬ 
ple,  passing  through  the  centerline  of  the  sample  in 
a  direction  parallel  to  the  ambient  wind  (see  Fig.  1). 


SAMPLE  WIDTH  EFFECTS  ON  MICROGRAVITY  FLAME  SPREAD 
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Thus,  in  the  three-dimensional  calculations,  one- 
quarter  of  a  corresponding  experimental  volume  was 
simulated.  The  physical  dimensions  of  the  three-di¬ 
mensional  computational  domain  were  12:8  X  6.4 
X  4.8  cm  in  the  x,  y,  and  z  directions,  respectively 
(12.8  X  4.8  cm  for  two-dimensions);  the  overall 
sample  size  in  the  three-dimensional  simulations  was 
no  larger  than  10  cm  long  X  4.5  cm  wide.  Note  that 
taking  into  account  the  symmetry  assumptions,  this 
corresponds  to  simulating  flame  spread  on  a  10  cm 
long  X  9  cm  wide  sample. 

The  numerical  grid  was  uniform  in  the  horizontal 
directions  and  stretched  in  the  vertical.  Cell  sizes 
were  0.66  mm  in  the  horizontal  directions  (parallel 
to  the  sample);  they  start  at  0.25  mm  at  the  sample 
surface  and  stretch  to  2  mm  at  the  top  of  the  com¬ 
putational  domain.  Various  grid  resolutions  were 
used  to  ensure  grid  independence  of  the  results 
(e.g.,  flame  spread  rates,  quenching  conditions).  A 
typical  three-dimensional  grid  contained  about 
600,000  grid  points  (192  X  96  X  32).  For  the  three- 
dimensional  runs,  approximately  14  h  CPU  time  are 
needed  for  1  s  of  simulated  real  time  on  current- 
generation  workstations. 

The  radiative  flux  on  the  sample  ramped  up  in 
time  following  tanh[t/(0.75  s)];  at  t  =  3.5  s  it  was 
turned  off.  This  procedure  is  consistent  with  previ¬ 
ous  simulation  [3,6]  and  experimental  studies  [9], 
This  flux  decayed  exponentially  with  distance  from 
the  x  =  0  as  exp[  —  (x//i)2].  The  half- width  profile  h 
=  0.25  cm.  Convective  oxygen  supply  to  the  flame 
was  controlled  by  varying  the  magnitude  of  the  am¬ 
bient  wind  uw  which  was  composed  of  33%  oxygen 
(Y0  =  0.33).  With  Y0  =  0.33,  the  ambient  wind 
speed  corresponding  to  an  oxygen-limited  regime 
was  known  from  the  authors’  previous  work  [6], 

Zero  gradient  conditions  are  used  for  the  farfield 
boundary  values  of  the  mass  fractions  and  tempera¬ 
ture  which  decay  exponentially  with  distance  from 
the  flame.  This  is  not  appropriate  for  the  velocity 
field  which  decays  quadratically  [8],  Instead,  bound¬ 
ary  conditions  were  obtained  from  an  analytically  de¬ 
termined  velocity  potential.  Following  Ref.  [8],  the 
velocity  on  a  computational  boundary,  iq,c,  is 

«bc  =  He  +  «w  =  V0  +  Mw, 

=  q  (2) 

Here,  ue  =  (ue,  ve,  we)  is  the  expansion  velocity  vec¬ 
tor,  »w  =  (mw,  0,  0)  is  the  ambient  wind  velocity 
vector  (see  Fig.  1),  0  is  the  remainder  potential,  q 
is  the  volumetric  heat  release  rate  of  combustion, 
are  /?„  and  h„  are  the  ambient  density  and  enthalpy. 
The  solution  for  0,  over  the  semi-infinite  domain  (s 
>  0)  can  be  calculated  using  a  Green’s  function,  G(x, 
y,  z ;  x',  (/',  s'), 


0(x,  y,  s,  t)  = 


J  |  V20  Gdx'dy'dz' 

+  J  |  0 - -  0)  G  dx'dy' 

where  on  the  sample  surface 

d0(x,  !/,  0)  3T 

Pochoo -  —  h  rii]  —  k  — 

3s  3s 


(3) 

(4) 


Here,  lii"  is  the  total  mass  flux  from  the  sample,  and 
k  is  the  thermal  conductivity.  Thus,  0  and,  therefore 
«bc  depend  on  the  mass  flux  and  heat  transfer  at  the 
sample  surface  and  on  the  chemical  heat  released  by 
the  gas-phase  reactions.  The  domain  is  open  to  the 
ambient  environment,  which  results  in  a  constant 
thermodynamic  pressure.  Note  that  in  an  experi¬ 
mental  study,  enclosure  effects  may  be  present 
which  are  not  accounted  for  in  the  simulation. 


Results  and  Discussion 

The  effects  of  a  finite  sample  width  on  flame  be¬ 
havior  were  investigated  by  simulating  samples  of 
three  different  widths:  2  cm,  4  cm,  and  9  cm.  To 
ensure  that  there  were  no  confounding  factors  re¬ 
lated  to  numerical  boundary  conditions,  the  size  of 
the  computational  domain  was  kept  the  same  for  sill 
cases.  In  experimental  studies  [9],  the  cellulosic  sam¬ 
ple  is  held  in  place  by  a  0.5  mm  thick,  American 
Iron  and  Steel  Institute  (AISI)  304,  stainless  steel 
holder.  Most  of  the  three-dimensional  simulations, 
therefore,  had  the  cellulosic  sample  bordered  by  a 
solid  with  the  thermal-physical  properties  of  AISI 
304  stainless  steel  (see  Fig.  1).  In  order  to  test  the 
effects  of  heat  loss  to  a  steel  sample  holder,  addi¬ 
tional  simulations  were  run  with  a  sample  holder 
which  had  the  thermal-physical  properties  of  the 
sample  (but  was  inert). 

This  paper  focuses  on  the  effects  of  a  finite  width 
sample  on  flame  behavior  during  two  different  re¬ 
gimes:  the  transition  period  between  localized  igni¬ 
tion  and  flame  spread  and  the  posttransition  period 
of  established,  possibly  steady-state,  flame  spread. 
To  survive  transition,  a  sufficient  amount  of  oxygen 
must  reach  the  reaction  zone.  In  these  simulations 
(Yq  =  0.33),  the  oxygen  supplied  by  an  ambient 
wind  of  uw  1.5  cm/s  or  below  was  insufficient  for 
the  flame  to  survive  transition  [6],  Three  different 
ambient  wind  conditions  were  used,  uw  =  2  cm/s, 
5  cm/s,  10  cm/s.  When  uw  =  2  cm/s,  flame  transi¬ 
tion  and  spread  occurs  in  an  oxygen-limited  environ¬ 
ment.  In  this  case,  the  flame  is  expected  to  be  rela¬ 
tively  weak  and,  therefore,  more  sensitive  to  changes 
in  its  environment  (such  as  heat  losses  to  a  steel  ver¬ 
sus  a  paper  sample  holder).  As  uw  increases  from  2 
cm/s,  the  oxygen  supply  improves  and  the  flame  be¬ 
comes  more  robust.  At  «w  values  larger  than  those 
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Fig.  2.  Flame  location  along  the  centerline  versus  time 
for  all  three-dimensional  cases  with  a  stainless  steel  holder 
(no  cases  with  a  paper  holder  are  shown).  Two  two-dimen¬ 
sional  cases  are  also  shown,  both  as  dotted  lines.  The  step¬ 
like  behavior  seen  on  the  curves  is  not  physical,  but  a  con¬ 
sequence  of  solving  the  discretized  equations. 

used  here,  heat  loss  to  the  ambient  gas  will  eventu¬ 
ally  extinguish  the  flame  (blowoff)  [1], 

Effects  of  Finite  Sample  Width  on  Transition  to 
Flame  Spread 

In  Fig.  2,  the  location  of  the  centerline  flame 
front,  as  determined  by  the  Yc  =  0.1  level  of  char 
mass  fraction,  is  plotted  versus  time  as  the  flame 
spreads  upwind.  The  trend  seen  in  Fig.  2  is  that  as 
the  ambient  wind  speed  increased  from  itw  =  2  cm/ 
s  to  the  maximum  of  «w  =  10  cm/s,  the  flame  spread 
rate  increased.  This  is  consistent  with  earlier  exper¬ 
imental  and  numerical  findings  [1,3,6],  Also,  as  the 
sample  width  increased  from  d  =  2  cm  to  a  maxi¬ 
mum  of  d  —  9  cm,  the  spread  rate  decreased.  At 
early  times  a  clear  period  of  transition  (non¬ 
monotonic  flame  spread)  existed  only  for  the  ww  = 
2  cm/s  cases.  The  duration  of  this  transition  period, 
from  longest  to  shortest,  was  approximately  5  s  (two- 
dimensional  case),  4  s  (three-dimensional,  d  =  4  cm 
and  9  cm  cases),  and  3.4  s  (three-dimensional, 
d  =  2  cm,  case);  note  that  radiative  flux  was  turned 
off  at  t  =  3.5  s.  Thus,  the  duration  of  the  transi¬ 
tion  period  increased  with  sample  width  when 
mw  =  2  cm/s. 

As  discussed  above,  the  ttw  =  2  cm/s  cases  were 
in  an  oxygen-limited  environment.  This  implies  that 
flame  transition  will  be  relatively  sensitive  to  changes 
in  the  net  inflow  velocity.  Thermal  expansion  from 


Fig.  3.  Top  view  of  flame  shape,  convective  (on  bottom 
half  of  each  figure)  and  net  (on  top  half  of  each  figure) 
oxygen  mass  flux  vectors,  and  color  contours  of  the  sample 
temperature  are  shown  for  the  »w  =  2  cm/s,  2  cm  wide 
sample,  (a)  Paper  sample  holder,  (b)  Steel  sample  holder. 
The  mass  flux  vectors  are  in  a  horizontal  plane  at  z  =  2.5 
mm  above  the  sample  surface.  The  time  is  t  =  20  s.  Inflow 
is  from  right  to  left.  The  entire  width  and  a  portion  of  the 
sample  length  is  shown  with  1  cm  spaced  grid  lines.  The 
flame  shape  is  presented  as  the  50  /ig/(cm3  s)  reaction  rate 
isosurface  (in  red).  The  oxygen  mass  flux  vectors  represent 
the  convective  only  (no  diffusion)  mass  flux  of  oxygen.  This 
quantity  was  used  in  order  to  mimic,  approximately,  use  of 
upstream  smoke  injection  for  flow  visualization,  as  was 
done  in  the  experiments  of  Ref.  [10]. 

the  flame  reduced  the  ambient  wind  by  an  amount 
(«e  <  0)  determined  by  the  analytical  solution  for 
the  potential  flow  (see  equation  2).  The  net  inflow 
velocity,  un  —  uw  +  ue,  was  used  in  the  simulations 
as  the  inflow  boundary  condition  in  the  x-direction 
(i.e.,  un  is  the  x  component  of  iq, c  in  equation  2). 

The  effects  of  thermal  expansion  on  the  incoming 
flow  have  been  visualized  experimentally  using  laser- 
illuminated  smoke  filaments  in  flame  spread  over  a 
liquid  pool  in  microgravity  [10],  The  convective  ox¬ 
ygen  mass  flux  vectors  on  the  bottom  half  of  Fig.  3 
show  that  the  incoming  flow  circumvented  the  flame 
in  a  manner  very  similar  to  that  seen  in  Ref.  [10], 
The  incoming  flow  speed  toward  the  flame  front  was 
dramatically  reduced  (note  size  of  arrowheads)  by 
thermal  expansion  in  both  the  paper  sample  holder 
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Fig.  4.  Top  view  of  both  the  9  cm  (a)  wide  and  the  2  cm 
wide  (b)  samples,  with  uw  =  2  cm/s  and  paper  sample 
holder.  The  time  is  t  =  10  s.  Streamlines  are  shown.  The 
50  /ig/(cm3  s)  reaction  rate  isosurface  is  color  contoured 
with  temperature  values.  The  grid  in  the  background  is 
composed  of  I  cm  spaced  lines  and  does  not  correspond 
to  the  sample  boundaiy  except  along  the  upwind  edge. 


case  (Fig.  3a)  and  the  steel  sample  holder  case  (Fig. 
3b).  Fig.  4  shows  the  flowfield  environment  for  the 
9  cm  wide  (Fig.  4a)  and  the  2  cm  wide  (Fig.  4b) 
sample  cases.  It  is  clear  that  in  the  2  cm  wide  sample 
case,  the  flame  was  immersed  in  a  three-dimensional 
flowfield,  even  in  the  centerline  region  since  the 
streamlines  diverged  laterally.  In  the  9  cm  wide  case, 
there  was  very  little  lateral  divergence  of  the  stream¬ 
lines  near  the  centerline,  suggesting  that  in  the  cen¬ 
terline  area  the  flow  is  nearly  two-dimensional.  The 
temperatures  of  the  simulated  flame  seen  in  Fig.  4 
are  in  the  range  measured  in  experiments,  1250  K 
to  1450  K  (thermocouple  located  2  mm  above  the 
sample  surface)  [11], 

Figure  5  is  a  plot  of  the  net  inflow  velocity,  un,  at 
the  centerline  versus  time  for  the  uw  =  2  cm/s, 
three-dimensional  cases  with  widths  d  —  2  cm,  4 
cm,  and  9  cm;  the  two-dimensional  case  is  also 
shown.  The  volumetric  heat  release  rate  of  combus¬ 
tion  (q),  which  was  used  in  the  computation  of  un, 
was  integrated  over  the  flame  volume  and  plotted 
versus  time  for  the  three-dimensional  cases  in  Fig. 
6.  Here,  q  was  determined  by  summing  the  product 
[heat  of  combustion]  X  [local  cell  volume]  X  [re¬ 
action  rate]  for  all  computational  cells  in  which  the 
reaction  rate  was  greater  or  equal  to  5  X  10  ~ 5 


Fig.  5.  Net  inflow  velocity  versus  time  for  the  three- 
dimensional  simulations  of  widths  d  =  2  cm,  4  cm,  and  9 
cm  and  ambient  wind  u„  =  2  cm/s.  The  two-dimensional 
uw  =  2  cm/s  case  is  also  shown.  These  values  of  the  inflow 
velocity  are  at  a  specific  height  (z  =  2.5  mm)  in  the  free- 
stream  (upwind  of  the  sample)  above  the  centerline  in  the 
inflow  boundary  plane  of  the  simulation.  They  are  char¬ 
acteristic  of  other  vertical  locations  above  the  centerline. 
Rapid,  ignition-induced,  thermal  expansion  at  early  times 
can  be  clearly  seen  in  all  cases. 

g/(cm3  s).  Rapid,  ignition-induced  heat  release  be¬ 
gan  at  f  =“  1  s  and  caused  a  dramatic  decrease  in  the 
net  inflow  velocity.  This  decrease  was  largest  for  the 
widest  flame  (the  two-dimensional  case)  and  small¬ 
est  for  the  narrowest  flame  (2  cm  wide  three-dimen¬ 
sional  flame).  After  a  sufficient  time  (the  transition 
period),  all  the  flames  reached  a  relatively  steady 
state.  Note  that  the  duration  of  the  transition  period 
(in  both  Fig.  5  and  Fig.  2)  scales  with  the  how  much 
the  ignition  event  decreased  the  inflow  velocity  rela¬ 
tive  to  its  ambient  value,  uw. 

Figure  6  also  shows  the  heat  release  rate  for  the 
three-dimensional,  9  cm  wide  flames  at  the  higher 
ambient  wind  velocities,  uw  =  5  cm/s  and  10  cm/s. 
The  ignition-induced  peak  heat  release  rate  was 
nearly  identical  for  all  three  9  cm  wide  cases.  For 
this  reason,  the  ignition-induced  expansion  velocity 
was  also  nearly  identical,  ue  =  —  0.4  cm/s.  This  was 
also  true  for  the  two-dimensional  cases  with  differ¬ 
ent  uw  (not  shown),  ue  “  —  1  cm/s.  Thus,  for  a  given 
flame  width,  the  ignition-induced  ue  is  independent 
of  the  ambient  wind  velocity,  uw.  This  explains  why, 
as  uw  increases,  the  duration  of  transition  decreases. 
Not  only  does  the  oxygen  supply  improve  as  «w  in¬ 
creases  (making  for  a  more  robust  flame),  but  the 
initial,  ignition-induced,  expansion-induced  velocity 
is  a  relatively  smaller  perturbation  from  the  ambient 
wind. 
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Fig.  6.  Heat  release  rate  of  gas-phase  reaction  versus 
time  for  all  the  three-dimensional  cases  with  a  stainless 
steel  holder.  The  heat  release  rate  was  found  by  summing 
q  (see  equation  2)  over  all  computational  cells  for  which 
the  gas-phase  reaction  rate  was  1  fig! (cm3  s)  or  greater. 
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Fig.  7.  Upstream  centerline  flame  spread  rate  versus 
width,  d,  of  the  sample  for  an  ambient  wind  of  uw  =  2  cm/ 
s.  Symbols  correspond  to  three-dimensional  cases.  The 
case  for  d  =  4  cm  with  a  paper  sample  holder  is  not  shown. 
Flame  spread  rate  for  two-dimensional  case  is  shown  as  a 
dashed  line. 
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There  was  essentially  no  difference  in  the  transi¬ 
tion  behavior  of  the  steel  and  paper  sample  holders 
(ww  =  2  cm/s  and  d  —  2  cm  and  9  cm).  Differences 
could  have  arisen  due  to  greater  heat  loss  to  the  steel 
sample  holder.  Apparently,  these  effects  were  not 


significant  during  transition.  They  did,  however,  in¬ 
fluence  the  post-transition  spread  rate  for  the  d  = 
2  cm  case,  as  shown  below. 

These  results  show  that  the  sample  size  has  an 
effect  on  the  duration  of  the  transition  period,  es¬ 
pecially  at  low  ambient  wind  speeds.  Some  cases 
have  transition  periods  that  are  comparable  to  the 
duration  of  reduced  gravity  in  drop  towers.  While 
the  results  of  these  simulations  have  not  yet  been 
validated  against  experimental  results,  they  suggest 
that  choosing  a  microgravity  facility  should  be  done 
with  these  issues  in  mind. 

Effects  of  Finite  Sample  Width  on  Established 
Flame  Spread 

Flame  spread  rates  were  used  to  investigate  the 
effects  of  a  finite  sample  width  on  established  flame 
spread.  These  were  computed  by  taking  the  slope  of 
the  flame  spread  curve  (e.g.,  in  Fig.  2)  over  a  time 
interval  in  which  the  flame  was  developing  in  a  rela¬ 
tively  steady-state  manner.  Some  care  was  needed  in 
determining  this  time  interval.  The  minimum  time 
used  was,  fmin  =  5  s,  which  was  after  the  transition 
period  for  all  flames  (e.g..  Fig.  2).  The  maximum 
time  was  chosen  to  be  before  the  presence  of  any 
confounding  effects  due  to  proximity  to  the  numer¬ 
ical  boundary.  Only  the  fastest  traveling  flames  (the 
three-dimensional  and  two-dimensional  cases  with 
«w  =  5  cm/s  and  10  cm/s)  were  affected  by  the 
numerical  boundary  (for  times  greater  than  approx¬ 
imately  15  s). 

Figure  7  shows  the  flame  spread  rates,  with  uw  = 
2  cm/s,  versus  width  for  the  three-dimensional  cases. 
Results  for  both  the  steel  and  paper  sample  holders 
are  plotted.  The  spread  rate  for  the  two-dimensional 
case  is  shown  as  a  dashed  line.  This  figure  shows  that 
as  the  sample  width  increased,  the  flame  spread  rate 
for  both  sample  holders  approached,  but  did  not 
reach,  that  of  the  two-dimensional  flame.  When  d 
=  2  cm,  the  flame  spread  rate  was  significantly 
slower  for  the  case  of  a  steel  versus  a  paper  sample 
holder.  This  implies  that  relatively  large  heat  losses 
to  the  steel  sample  holder  slowed  the  flame  since 
the  spread  rates  were  nearly  identical  for  the  widest 
(d  =  9  cm)  sample  case.  This  interpretation  is  cor¬ 
roborated  by  the  fact  that  flame  spread  in  the  d  = 
4  cm  case  is  as  fast  as  the  d  =  2  cm  case,  even 
though,  based  on  the  net  inflow  velocity  in  Fig.  5, 
the  d  —  2  cm  case  had  a  better  oxygen  supply.  Con¬ 
sider  also  Fig.  3,  which  shows  the  temperature  of 
the  sample  for  both  the  paper  (Fig.  3a)  and  steel 
(Fig.  3b)  sample  holder  cases.  Heat  loss  effects  are 
clearly  greater  in  Fig.  3b:  the  flame  was  smaller,  the 
sample  temperature  distribution  was  more  confined 
to  the  center  of  the  sample,  and  the  peak  sample 
temperature  was  lower.  The  distinct  difference  in 
the  temperature  distribution  between  Fig.  3a  and  b 
was  due  to  a  relatively  strong  (hot)  flame  above  the 
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Fig.  8.  Upstream  centerline  flame  spread  rate  versus  the 
ambient  wind  speed,  »w.  Three-dimensional  cases  with  the 
steel  sample  holder  are  shown  along  with  two-dimensional 
cases. 

sample  edge  in  Fig.  3a.  The  net  oxygen  mass  flux 
vectors  on  the  upper  half  of  the  figures  show  that 
this  edge  flame  consumed  the  oxygen  diffusing  from 
the  outer  free  stream  before  it  reached  the  center. 
No  such  edge  flame  could  exist  in  steel  sample  case 
since  heat  losses  at  the  sample/steel  boundary  pre¬ 
vented  the  sample  from  reaching  the  pyrolysis  tem¬ 
perature.  Note  that  the  point  of  Fig.  7  is  that  even 
with  a  sample  sufficiently  wide  to  ensure  that  edge 
effects  were  negligible  on  the  spread  rate  of  the  cen¬ 
ter  flame,  the  spread  rate  was  still  not  that  of  the 
two-dimensional  flame. 

Figure  8  is  a  plot  of  spread  rate  versus  ambient 
wind  speed  from  the  three-dimensional  9  cm  wide, 
steel  holder  cases  and  the  two-dimensional  cases.  It 
is  clear  that  as  the  ambient  wind  speed  increased, 
the  spread  rate  of  the  finite  width,  9  cm  wide  flame 
tended  to  the  two-dimensional  value.  This  occurred 
for  the  same  reason  the  duration  of  transition  de¬ 
creased  as  »w  increased:  with  increased  oxygen  sup¬ 
ply  a  given  width  flame  was  less  influenced  by  other 
environment  factors.  For  example,  the  net  inflow  ve¬ 
locities  in  the  9  cm  three-dimensional  and  the  two- 
dimensional  cases  differed  by  0.5  cm/s  (not  shown) 
when  uw  =  10  cm/s.  This  is  5  times  their  difference 
when  uw  =  2  cm/s  (see  Fig.  5).  Bnt  three-  and  two- 
dimensional  flame  spreads  agreed  when  »w  =  10 
cm/s  because  there  was  an  abundance  of  oxygen. 
This  can  easily  be  seen  in  Fig.  8  by  the  relative  in¬ 
sensitivity  of  the  flame  spread  rate  to  changes  in  uw 
(i.e.,  the  formation  of  a  plateau)  as  uw  increases.  Al¬ 
though  the  simulations  were  not  performed,  it  is 
possible  that  even  for  the  2  cm  sample  width,  the 


three-  and  two-dimensional  flame  spread  rates 
would  be  in  better  agreement  as  u„  increased. 

Conclusions 

The  width  of  the  sample  affected  both  the  dura¬ 
tion  of  the  flame  transition  period  and  the  post-tran¬ 
sition  flame  spread  rate.  Finite  width  affects  were 
most  significant  when  the  ambient  wind  was  rela¬ 
tively  small  (limited  oxygen  supply).  In  such  envi¬ 
ronments,  the  velocity  due  to  thermal  expansion 
from  the  flame  reduced  the  net  inflow  of  oxygen 
enough  to  significantly  affect  flame  behavior.  The  ex¬ 
pansion  velocity  increased  with  increasing  flame  size 
(i.e.,  sample  width).  For  a  given  sample  width,  the 
influence  of  thermal  expansion  on  the  net  incoming 
oxygen  supply  decreased  as  the  ambient  wind  speed 
increased.  Thus,  both  the  transition  and  flame 
spread  behavior  of  the  three-dimensional  flame 
(along  the  centerline)  tended  to  that  of  the  two-di¬ 
mensional  flame  with  increasing  ambient  wind 
speed.  These  results  emphasize  the  importance  of 
including  velocity  boundary  conditions  which  ac¬ 
count  for  thermal  expansion  effects  at  low  ambient 
wind  velocities. 

Heat  losses  to  the  sample  holder  were  found  to 
significantly  lower  the  temperature  on  the  sample 
near  the  sample-sample  holder  interface.  As  a  re¬ 
sult,  the  flame  spread  rate  in  the  case  of  the  narrow¬ 
est  sample,  with  the  slowest  ambient  wind,  was  re¬ 
duced.  These  results  are  based  on  simulations  only. 
Further  microgravity  experiments  are  needed  to 
conclusively  validate  these  findings. 
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